Topological insulators with SU(2) Landau levels 
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We construct continuum models of 3d and 4d topological insulators by coupling spin-i fermions 
to an SU(2) background gauge field, which is equivalent to a spatially dependent spin-orbit coupling. 
Higher dimensional generalizations of flat Landau levels are obtained in the Landau-like gauge. The 
2d helical Dirac modes with opposite helicities and 3d Weyl modes with opposite chiralities are 
spatially separated along the third and fourth dimensions, respectively. Stable 2d helical Fermi 
surfaces and 3d chiral Fermi surfaces appear on open boundaries, respectively. The charge pumping 
in 4d Landau level systems shows quantized 4d quantum Hall effect. 
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Time-reversal (TR) invariant topological insulators 
(TI) have become a major research focus in condensed 
matter physics[l~3]. Different from the 2d quantum Hall 
(QH) and quantum anomalous Hall systems which are 
topologically characterized by the first Chern number [jj- 
6], time- reversal invariant TIs are characterized by the 
second Chern number in 4d 0, H| and the index in 2d 
and 3d [Uli | • Various 2d and 3d TIs are predicted the- 
oretically and identified experimentally exhibiting stable 
gapless ID helical edge and 2d surface modes against TR 
invariant perturbations 11, T^-l^. Topological states 



have also been extended to systems with particle-hole 
symmetry and superconductors 2(11- 22j|. 

Most current studies of 2d and 3d TIs focuses on Bloch- 
wave bands in lattice systems. The non-trivial band 
topology arises from spin-orbit (SO) coupling induced 
band inversions 11] . However, Landau levels (LL) are es- 
sential in the study of QH effects because their elegant 
analytical properties enable construction of fractional 
states. Generalizing LLs to high dimensions gives rise to 
TIs with explicit wavefunctions in the continuum, which 
would facilitate the study of the exotic fractional TIs. Ef- 
forts along this line were pioneered by one of the authors 
and Hu |7|. They constructed LLs on the compact S 4 
sphere by coupling fermions with the SU(2) monopole 
gauge potential, and various further developments ap- 
peared [23l - [27| . 2d TIs based on TR invariant LLs have 
also been investigated [l(J. Two of the authors general- 
ized LLs of Schrodinger fermions to high-dimensional flat 
space |28] by combining the isotropic harmonic potential 
and SO coupling. LLs have also been generalized to high 
dimensional Dirac fermions and parity-breaking systems 

urn 

In all the above works, angular momentum is explic- 
itly conserved, thus they can be considered as LLs in 
the symmetric-like gauge. In 2d, LL wavefunctions in 
the Landau gauge are particularly intuitive: they are 
Id chiral plane- wave modes spatially separated along the 
transverse direction. The QH effect is just the Id chi- 



ral anomaly in which the chiral current generated by the 
electric field becomes the transverse charge current. In 
this article, we develop high dimensional LLs with flat 
spectra as spatially separated helical Dirac or chiral Weyl 
fermion modes, i.e., the SU(2) Landau-like gauge. They 
are 3d and 4d TIs defined in the continuum possessing 
stable gapless boundary modes. For the 4d case, they 
exhibit the 4d quantum Hall effect which is a quan- 
tized non-linear electromagnetic response related to the 
spatially separated (3+l)d chiral anomaly. Our methods 
can be easily generalized to arbitrary dimensions and also 
to Dirac fermions. 

The 3d case We begin with the 3d TR invariant LL 
Hamiltonian for a spin-^ fermion as 



H 



3d,i/=± 
LL 



2m 



-muj 2 so z 2 



VLU so Z(p x (Jy -Py<T X )(l) 



which couples the ID harmonic potential in the z di- 
rection and the 2d Rashba SO coupling through a z- 
dependent SO coupling strength. H^£ v possess trans- 
lation symmetry in the xy-plane, TR and parity sym- 
metries. H^ L ,+ and H^ d £~ share the same physics up 
to a z — > —z transformation. For simplicity, unless 
explicitly mentioned, we only consider H^ d £ + and sup- 
press the index v. Eq. [1] can be reformulated in the 
form of an SU(2) background gauge potential as Hjf L = 
^(p- - C A) 2 - \muj 2 z 2 , where u> so = \eG\/mc, and A 
takes the Landau-like g aug e as A x = Ga y z, A y = —Ga x z 
and A z = 0. In Ref. [28], a symmetric-like gauge with 
A' = Go x r is used, which explicitly preserves the 3d 
rotational symmetry. However, the SU(2) vector poten- 
tials A' and A are not gauge equivalent. As shown below, 
the physical quantities of Eq. [I] such as density of states 
(DOS), are not 3d rotationally symmetric. Nevertheless, 
these two Hamiltonians belong to the same topological 
class. A related Hamiltonian is also employed for study- 
ing electromagnetic prop erties in superconductors with 
cylindrical geometry [3l| . 

Eq. ([1]) can be decomposed into a set of ID harmonic 
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oscillators along the z-axis exhibiting flat spectra, a key 
feature of LLs. We define a characteristic SO length scale 

— . Each of the reduced ID harmonic oscil- 



^so 

hit or Hamiltonian is associated with a 2d helical plane- 
wave state as H z (k 2d ) = ^+\muj 2 s J [ z~l 2 so k 2d 2Z 2d (k 2d )\ 2 , 
where k 2d — (k 2 + k 2 )? and k 2 d = {k x ,k y ); the helicity 
operator is defined as T< 2 d(k 2 d 
LL eigenstates are solved as 



k x Oy 



k y <7 x . The n-th 



J ^n[«-«o(A;2d,S)]®XE(A;2d), (2) 



where r 2d — (x, y); XT,{k 2 d) are eigenstates of the helicity 
satisfying z^XT,{k 2 d) — £xs(^2d) with helicity eigenvalues 
E = ±1; (f> n [z — Zq] are the eigenstates of the n-th har- 
monic levels with the central positions located at zo, and 
zo{k 2 d, E) = Yil 2 so k 2 d- The energy spectra of the n-th LL 
is E n — (n + i)fiw so , independent of k 2 d and E. 
The classical equations of motion for H 

as 



£^ are derived 



r = — , p 2d = 0, p z = 2to(p x S) z - muj 2 z, 



rn 
2lj 



">2</ 



-zS z p 2d , 



S z 



2lj 



-zS 2d -p 2 d, 



(3) 



where S^d = (gOx, \(Jy) and S 2 = A<7z. If we choose 
the initial condition of spin S such that S Zt o = and 
S 2 d : o JL P2d,0; then 5 is conserved and lies in the xy-plane. 
The motion is reduced to a coplanar cyclotron one in the 
vertical plane perpendicular to S. In other words, the 
orbital angular momentum and spin are locked, a feature 
from SO coupling. 
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FIG. 1: (Color online) 3d and 4d LLs for and Kf L as 

spatially separated 2d helical SO plane-wave modes localized 
along the z-axis (A), and 3d Weyl modes localized along the u- 
axis (B), respectively. Their central locations are zo(fc2d, S) = 
S/sofed and wo(fe3d,S) = T,ls k3d, respectively. Note that 2d 
Dirac modes with opposite helicities and the 3d ones with 
opposite chiralities are located at opposite sides of z = and 
u — planes, respectively. 

In the 2d LL case, the spatial coordinates x and y 
are non-commutative if projected to a given LL, say, 
the lowest LL (LLL). The LLL wavefunctions are chi- 
ral Id plane-waves along the x-direction whose central 



y-position is yo = l%k x , thus Id chiral modes with oppo- 
site chiralities are spatially separated. Consequently, the 
xy-plane can be viewed as the 2D phase space of a Id 
system, in which y plays the role of k x . The momentum 
cut-off of the bulk states is determined by the system size 
along the y-direction as IfcJ < Similarly, the 3d LL 
wavefunctions in Eq. [5] are spatially separated 2d helical 
plane- waves along the z-axis. As shown in Fig. [1] (A), 
for states with opposite helicity eigenvalues, their central 
positions are shifted in opposite directions. After the pro- 
jection into the LLL, z is equivalent to l 2 a (k x a y — k y a x ), 
and thus 

[x, z] LLL = il 2 so (T y , [y, z] LLL = -il 2 so <T x , [x, y] LLL = 0. (4) 

Interestingly, the 3d LL states can be viewed as states 
in the 4D phase space of a 2d system (with x and y 
coordinates), in which |z| plays the role of the magnitude 
of the 2d momentum. In fact, the momentum cut-off of 
the bulk states is exactly determined by the spatial size 
L z along the z direction as 



k 



2d 



<k B = 
K 2d — 



2l 2 



(5) 



in which -fy- < z < 1 



2 ^ *, ^ 2 . Applying the open boundary 
condition along the z direction, and periodic boundary 
conditions in the x and y direction with spatial sizes L x 
and L y respectively, we can easily count the total number 



of states to be Af = 



The L 2 dependence of N 



may seem puzzling for a 3d system, but expressing L z in 
terms of Eq. ([5]), we find N — ■^L x L y (kf d ) 2 , which is 
the conventional state counting of a 2d system expressed 
in terms of the 4D phase space volume. 

The topological property of this 3d LL system man- 
ifests through its helical surface spectra. Let us con- 
sider the upper boundary located at z = zb for H^ L . 
For simplicity, we only consider the LLL as an exam- 
ple. If zb > 0, the positive helicity states 'I'q fc 2d e=i 
with k 2 d > k^d = zbIJ 2 are confined at the bound- 
ary. The Id harmonic potential associated with k 2 d and 
E = 1 is truncated at z = zb, and thus the surface 
spectra acquire dispersion. If we neglect the zero-point 
energy, the surface mode dispersion is approximated by 
^muj 2 (k — kfd) 2 ^ with k > k® d . If the chemical poten- 
tial ^ lies above the LLL, it cuts the spectra at surface 
states with a Fermi wavevector kf > k^ d - The Fermi 
velocity is Vf ~ m(kf — k^^so^so- The surface Hamilto- 
nian is approximated as H s f s=s vj(p x a) ■ z — fi with an 
electron- like Fermi surface with E = 1. In another case, if 
the upper boundary is located at zb < 0, then the nega- 
tive helicity states *f? g 2 s = _i with k 2 d < k^ d , and all the 
positive helicity states are pushed to the boundary as sur- 
face modes. Depending on the value of p, we can have a 
hole-like Fermi surface with E = — 1, a Dirac Fermi point, 
or, an electron-like Fermi surface with E — 1. Similarly, 
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any other LL gives rise to a branch of gapless helical 
surface modes, and each filled bulk LL contributes one 
helical Fermi surface. For the lower boundary, the anal- 
ysis is parallel to the above. Each filled LL gives rise to 
an electron- like helical Fermi surface with E = — 1, or 
hole-like with E = 1. 

The 4d case The above procedure can be straightfor- 
wardly generalized to any higher dimension. For exam- 
ple, the 4d LL Hamiltonian is denoted as 

Tr4d,u=± Pu i 2 2, Pzd - - fc\ 

H L £ = h -mu) u + - vuup 3d ■ a, 6) 

2m 2 2m 

where u and p u refer to the coordinate and momen- 
tum of the 4th dimension, respectively, and p 3d is the 
3-momentum in the xyz-space. Without lose of the gen- 
erality, let us focus on the case of v = + and Eq. [5] can 
be represented as Hjf L = ^ Yn=i{Pi _ 7 A i) 2 ~ muj 2 u 2 , 
where the SU(2) vector potential takes the Landau-like 
gauge with Ai = GoiU for i = x,y,z and A u — 0. Eq. [6] 
preserves the translational and rotational symmetries in 
the xyz-sp&ce and TR symmetry. Similar to the 3d case, 
the 4d LL spectra and wavefunctions are solved by reduc- 
ing Eq. [6] into a set of Id harmonic oscillators along the 
w-axis as H u (k 3d ) = f£ + ^moj 2 (u - l 2 so k 3d i: 3d ) 2 where 
k 3d = {k 2 + k 2 + k 2 )? and t, 3d — fed • <?. The LL wave- 
functions are 

*„,fc 3d ,s( f ' w ) = j %Sd *<i>n[u- wo(fc 3 d,E)] ® Xe(M, (7) 

where, the central position M (fc 3( j,E) = Yil 2 so k 3d ; xs are 
eigenstates of 3d helicity E 3( ; with eigenvalues E = ±1. 
Inside each LL, the spectra are flat with respect to k 3d 
and E. This realizes the spatial separation of the 3d 
Weyl fermion modes as shown in Fig. [1] (B). After the 
LLL projection for Hjf L , the non-commutative relations 
among coordinates are [rt, u] — il 2 so cn and [n, rj] = for 
i = x,y, z. Similar to the discussions of the 3d case, the 
4d LLs can be viewed as states of the 6D phase space of 
a 3d system: increasing the width along the u-direction 
corresponds to increasing the bulk momentum cut-off 
k 3d < kf d = L u /(2l 2 so ). If the LLL is fully-filled, the 

total number of states is given by J\f = -tArf - " ■ Re- 
expressing L u = 2k 3d l 2 so we find Af = ^2L x L y L z (k^ d ) 3 , 
which is the conventional state counting of a 3d system 
expressed in terms of the 6D phase space volume. With 
an open boundary imposed along the u direction, 3d he- 
lical Weyl fermion modes appear on the boundary. 

Now let us consider the generalized 4d QHE[7|, |8| as 
the nonlinear electromagnetic response of (4 + l)d LL 
system to the external electric and magnetic (EM) fields, 
with E || B in the xyz-space. Without loss of generality, 
we choose the EM fields as E = Ez and B = Bz, which 
are minimally coupled to the spin- 1/2 fermion, 

Hf L {E,B) = -ttv 2 u + l -mu 2 (u + l l 2 D ■ a)\ (8) 




FIG. 2: The central positions uo(m, k z , v) of the 4d LLs in the 
presence of the magnetic field B = Bz. The branch of m = 
runs across the entire it-axis, which gives rise to quantized 
charge transport along M-axis in the presence of E \\ B as 
indicated in Eq. [TOl 

where D = V — i^-A em . Here A em is the U(l) mag- 
netic vector potential in the Landau gauge with A em x = 



0, A errii y = Bx and A errltZ = —cEt. We define Ib = 




where eB > is assumed. 

The £?-ficld further reorganizes the chiral plane-wave 
states inside the n-th 4d LL into a series of 2d magnetic 
LLs in the xy-plane. For the moment, let us set E z — 0. 
The eigenvalues E n — (n + ^)hoj so remain the same as 
before without splitting, while the eigen- wavefunctions 
are changed. We introduce a magnetic LL index m in the 
xy-plane. For the case of m — 0, the eigen- wavefunctions 
are spin polarized as 

^ n ,mMky,k z ) = e lk ^+ %k * z 4> n [u - u {k z , m = 0)] 

<E) Xm=o[x - x (k y )}, (9) 

where xo = l%k y and xo = [</>o{x — xo),0] T is the zero 
mode channel of the operator —iD ■ a with the eigenvalue 
of Aq = k z . The central positions of the it-direction har- 
monic oscillators are Uo(k z ,m = 0) = l 2 so k z . For m > 1, 
the eigen-modes of —iD ■ a come in pairs as Xm,±[% — 
xo{k v )] = [otm,±<t>m{x - x ), /3 m <j>m-i{x - x )] T , where, 
coefficients a m ,± = lBk z ± \Jl\k 2 -I- 2m, j3 m = —iy/2m, 

and the eigenvalues are A m ,± = ±\Jk 2 z + 2ml~ 2 . The 
corresponding eigen- wavefunctions 

e lk v y+lk * z (j) n [u-u (k z ,m,±)}x m ,±[x- x (ky)}, whe re the 

central positions uo(k Zl m, ±) = ±l 2 so \J k 2 + 2m( B 2 . 

For the solutions of Eq. [S]with the same 4d LL index n, 
the 2d magnetic LL with index m = is singled out. The 
central positions of states in this branch is linear with k z , 
and thus run across the entire u-axis, while those of other 
branches with m > 1 only lie in one half of the space as 
shown in Fig. [21 After turning on E z , k z is accelerated 
with time as k z (t) — k z (0) + eE z t/h, and thus the central 
position uo(m — 0, k z ) moves along the it-axis. Only the 
m = branch of the magnetic LL states contribute to 
the charge pumping which results in an electric current 
along the u-direction. Within the time interval At, the 
number of states with each filled 4d LL passing a cross 
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LyLy eEAtL z 
e dN 



section perpendicular to the it-axis is N 

which results in the electric current density L [ L ^ ■ If 
the number of fully filled 4d LLs is n occ , the total current 
density along the it-axis is 



4ir 2 h 



E-B, 



(10) 



where a = e 2 /(he) is the fine-structure constant. 
This quantized non-linear electromagnetic response is in 
agreement with results from the effective theory Q as the 
4d generalization of the QH effect. If we impose open 
boundary conditions perpendicular to the it-direction, 
the above charge pump process corresponds to the chi- 
ral anomalies of Weyl fermions with opposite chiralities 
on the two 3d boundaries, respectively. Since they are 
spatially separated, the chiral current corresponds to the 
electric current along the u-direction. 

Although the DOS of our 4d LL systems is different 
from the 4d TIs in the lattice system with translation 
symmetry Q, their electromagnetic responses obey the 
same Eq. [TUJ It is because that only the spin-polarized 
to = branch of LLs, '5 ' n ,m=o(k y , k z ), is responsible for 
the charge pumping. For this branch, B^-field quantizes 
the motion in the rry-plane so that the x and it coor- 
dinates play the role of k y and k z , respectively. Con- 
sequently, the system can be viewed as the 4D phase 
space of coordinates y and z, and scales uniformly as 
L x L y L z L u with the conventional thermodynamic limit 
of a 4d system. 

High dimensional LLs of Dirac electrons The above 
procedure can be generalized to Dirac fermions as square 
root problems to Eq. (JTJ and Eq. |6|). For the 3d case, 
we have 



H 



3d 

LL.Di 



Iso^sc 

~7T 







a ■ p - 



a -p+ifi 




(11) 



whose square exhibits a diagonal block form with a super- 



symmetric structure as (H[f irac ) 2 / (fiw M ) 



.H 



3d, 



2 i x± LL 

are E± T: 



f The eigenvalues of Eq. (dJ) 

±^/nhuj so . The eigen-wavefunctions are 
constructed based on the eigenstates of 2f^ 1/=± as 

_,±*~ . , J T , where 



^±n,D lr ac(k2d,^) 



n.k- 



' n-l,fe 2(i ,-S J 
z (k 2 d,^)} ® Xs{hd) is the 

eigenstate of H^ a £~. The 0th LL states are Jackiw-Rebbi 
half-fermion modes with only the upper two components 
The 4d Hamiltonian can be constructed 



\/2 L ~n,fc 2d .i;" 
%2d ' F2d <j>n[z + 
3d, 



nonzero 
as 



Ml33| 



H 



Ad 

LL,Di 



Iso^sc 





o ■ P3d + ij-, 

■I 2 , 



o -P3d~ ^- 3 




(12) 



where a u = — (u + i-j^Pu) is the phonon annihi- 
lation operator in the it-direction. The eigen- values 



are still E± 7 



^i*n,Dirac(k3d 



and the eigenstates are 



= ±y/nhuj so , 

si = 4=r^ + 

Laughlin-type wavefunctions The construction of the 
Laughlin-type wavefunctions for interacting fermions is 
difficult for these SO coupled high dimensional LL sys- 
tems. Nevertheless, for the 4d case in the magnetic field, 
the LLL states with both n — and m = are spin- 
polarized, and their total DOS is finite as p^L m = = 
47r2 p ;ii . Even though they are degenerate with other 

LLL states with (n — 0, m ^ 0), they are favored by 
repulsive interactions if they are partially filled. The 
Laughlin wavefunction in the Landau gauge for the 2d 
LLs has been constructed in Ref. 
to the 4d case, and define w = e 
then LL states are represented as w h:c v hz up to a Gaus- 



34] 



We generalize it 

and v = e l L * , 



sian factor e 2!so 21 b . The Laughlin-type wavefunction 
can be constructed as 

^(wi,vi;...;w NwNy ,v N;isNy ) = (det[w i *v i ']) q , (13) 

where q is an odd number; det[it;^ x i'^ z ] represents the 
Slater determinant for the fully filled single particle states 
w h *v h * and < h x < N x - 1 and < h z < N z - 1. 
The study of the topological properties of Eq. [T3] will be 
deferred to a later publication. 

In conclusion, we have generalized LLs into 3d and 4d 
in the Landau-like gauge by coupling spatially dependent 
SO coupling with harmonic potentials. This method can 
be generalized to arbitrary dimensions by replacing the 
Pauli-matrices with the T-matrices in corresponding di- 
mensions. These high dimensional LLs exhibit spatial 
separation of helical or chiral fermion modes with oppo- 
site helicities, which give rises to gapless helical or chiral 
boundary modes. The 4d LLs give rise to the quantized 
non-linear electromagnetic responses as a spatially sepa- 
rated (3+l)d chiral anomaly. 
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